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1. INTRODUCTION 

The present paper is devoted to degenerate electron plasma behaviour 
research. Analysis of processes taking place in plasma under effect of ex- 
ternal electric field, plasma waves oscillations with various types of condi- 
tions of electron reflection from the boundary has important significance 
today in connection with problems of such intensively developing fields as 
microelectronics and nanotechnologies [1] - [6] . 

The concept of "plasma" appeared in the works of Tonks and Langmuir 
for the first time (see [7]-[9]), the concept of "plasma frequency" was 
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introduced in the same works and first questions of plasma oscillations 
were considered there. However, in these works equation for the electric 
field was considered separately from the kinetic equation. 

A. A. Vlasov [10] for the first time introduced the concept of "self- 
consistent electric field" and added the corresponding item to the kinetic 
equation. Now equations describing plasma behaviour consist of anchor 
system of equations of Maxwell and Boltzmann. The problem of electron 
plasma oscillations was considered by A. A. Vlasov [10] by means of solution 
of the kinetic equation which included self-consistent electric field. 

L.D. Landau [11] had supposed that outside of the half-space contain- 
ing degenerate plasma external electromagnetic field causing oscillations 
in plasma is situated. By this Landau has formulated a boundary condi- 
tion on the plasma boundary. After that the problem of plasma oscillation 
turns out to be formulated correctly as a boundary problem of mathemat- 
ical physics. 

In [11] L.D. Landau has solved analytically by Fourier series the prob- 
lem of collisionless plasma behaviour in a half-space, situated in external 
longitudinal (perpendicular to the surface) electric field, in conditions of 
specular reflection of electrons from the boundary. 

Further the problem of electron plasma oscillations was considered by 
many authors. Full analytical solution of the problem is given in the works 
[12] and [13]. 

This problem has important significance in the theory of plasma (see, 
for instance, [2], [14] and the references in these works, and also [15], [16]). 

The problem of plasma oscillations with diffuse boundary condition was 
considered in the works [17], [18] by method of integral transformations. In 
the works [19], [20] general asymptotic analysis of electric field behaviour 
at the large distance from the surface was carried out. In the work [19] 
particular significance of plasma behaviour analysis close to plasma reso- 
nance was shown. And in the same work [19] it was stated that plasma 
behaviour in this case for conditions of specular and diffuse electron scat- 
tering on the surface differs substantially. 
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In the works [22] general questions of this problem solvability were con- 
sidered, but diffuse boundary conditions were taken into account. In the 
work [22] structure of discrete spectrum in dependence of parameters of 
the problem was analyzed. The detailed analysis of the solution in general 
case in the works mentioned above hasn't been carried out considering the 
complex character of this solution. 

The present work is a continuation of electron plasma behaviour in 
external longitudinal alternating electric field research [22] - [26] . 

In the present paper the linearized problem of plasma oscillations in slab 
in external alternating electric field is solved analytically. Specular bound- 
ary conditions for electron reflection from the boundary are considered. In 
[24]-[26] diffuse boundary conditions were considered. 

The coefficients of continuous and discrete spectra of the problem are 
obtained in the present work, which allows us to derive expressions for 
electron distribution function at the boundary of conductive medium and 
electric field in explicit form. 

The present work is a continuation of our work [27], in which questions 
of half-space plasma waves reflection from the plane boundary bounding 
degenerate plasma were considered. 

Let us note, that questions of plasma oscillations are also considered in 
nonlinear statement (see, for instance, the work [28], [29]). 

2. PROBLEM STATEMENT 

Let degenerate plasma occupy a slab (particularly, thin fllms) —a<x< 

a. 

We take system of equations describing plasma behaviour. As a kinetic 
equation we take Boltzmann — Vlasov r-model kinetic equation: 

-^ + v-^ + eE-^= -^^^^ ' ^ ' ' . 1.1 

at or op T 

Here / = /(r, v, t) is the electron distribution function, e is the electron 
charge, p = mv is the momentum of an electron, m is the electron mass, 
r is the character time between two collisions, E = E(r,i) is the self- 
consistent electric fleld inside plasma, feq = feqiy^ t) is the local equilibrium 
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Fermi — Dirac distribution function, feq = Q{£,Fit, x) — £), where 0(x) is 
the function of Heaviside, 



Six) = 



1, X > 0, 
0, x<0, 



E,p{t^x) = ^mvp{t,x) is the disturbed kinetic energy of Fermi, £ = ^mv'^ 
is the kinetic energy of electron. 

Let us take the Maxweh equation for electric field 

divE(r,t) = 47rp(r,t). (1.2) 
Here p(r, t) is the charge density, 

p(r, t) = ej (/(r, v, t) - /o(v)) dQp, (1-3) 

where 

2d p 

Here /o is the undisturbed Fermi — Dirac electron distribution function, 

/o(£) = e(£f-£), 

h is the Planck's constant, u is the effective frequency of electron collisions, 
= 1/r, £f = ^mv^p is the undisturbed kinetic energy of Fermi, vp is the 
electron velocity at the Fermi surface, which is supposed to be spherical. 

We assume that external electric field outside the plasma is perpendi- 
cular to the plasma boundary and changes according to the following law: 

Then one can consider that self-consistent electric field E(r, t) inside 
plasma has one a;-component and changes only lengthwise the axis x: 

E = {E^(x,t),0,0}. 

Under this configuration the electric field is perpendicular to the boun- 
dary of plasma, which is situated in the plane x = 0. 

We will linearize the local equilibrium Fermi — Dirac distribution f^q 
in regard to the undisturbed distribution /o(£): 
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where S{x) is the delta - function of Dirac. 

We also linearize the electron distribution function / in terms of absolute 
Fermi — Dirac distribution /o(£): 

/ = /o(£) + /i(:^,v,t). (1.4) 

After the linearization of the equations (1.1)-(1.3) with the help of (1.4) 
we obtain the following system of equations: 

^+^.|^+^/i(:^, V, t) = (5(£i.-£) [eE,{x, t)^,+^[£f (x, t) -£^^]] , (1.5) 
From the law of preservation of number of particles 

I f,,du, = I fdU, 

we find: 

[£i.(x,t) - 8,f] J Si^F - £.)d'p = J hd'p. (1.7) 

From the equation (1.5) it is seen that we should search for the function 
fi in the form proportional to the delta - function: 

h = £,F6i£,F-£^)Hix,fi,t), fi = ^. (1.8) 

V 

The system of equations (1.5) and (1.6) with the help of (1.7) and (1.8) 
can be transformed to the following form: 

dH dH 

~dt ^ ^^''^'dx ^ ^' ^ 

evFfi 



/I 
H{x,ii',t)dii', 



Further we introduce dimensionless function 

e{x,t) = ^^Ex{x,t) 

UCf 
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and pass to dimensionless coordinate xi = x/l , where I = vpT is the mean 
free path of electrons, and we introduce dimensionless time ti = ut. We 
obtain the following system of equations: 

dH dH 1 

-wr-^f^a \- i^H{xi,ii,ti) = fie{xi,ti) + - H{xi, /i' ,ti)dii' , (1.9) 

oil oxi 2 7-1 

= ^ H{xi, // , ti)d^ . (1.10) 
Here cjp is the electron (Langmuir) frequency of plasma oscillations, 

^ m 

N is the numerical density (concentration), m is the electron mass. 

We used the following well-known relation for degenerate plasma for the 
derivation of the equations (1.9) and (1.10) 



Let A; to be a dimensional wave number, and let us introduce dimension- 
less wave number ki = k—, then kx = where e = —. We introduce 



CUp € CUp 



1 • ^ 

the quantity Ui = ur = —. 

V 

3. BOUNDARY CONDITIONS STATEMENT 

Let us outline the time variable of the functions H{xi, /i, ti) and e{xi, ti), 
assuming 

^f(xi,//,ti) =e-^'^^*^/i(xi,//), (2.1) 

e{xi,ti) = e-'^^'^'eixi). (2.2) 

The system of equations (1.9) and (1.10) in this case will be transformed 
to the following form: 

1 

X f 

'^'dx~ ^ ^"^ ~ i^i)h{xi, fi) = iie{xi) + 2 / ^(^i' (^-3) 

-1 
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-1 

Further instead of xi, ti we write rc, t. We rewrite the system of equations 
(2.3) and (2.4) in the form: 

1 

dh If 

'^dx ^ '^o^^^' ^ ^^^^^ ^ 2 / ^^"^^ 

-1 

1 

'^''^'^^ ^ ^ h{x,fi')dfi'. (2.6) 



Here 



(ix 2^2 7 
-1 



Zq = 1 — iuji = 1 = 1 — iujT = 1 — ir2, = ujT. 

V 



We consider the external electric field outside the plasma is perpendi- 
cular to the plasma boundary and changes according to the following law: 
£^oexp(— ia;t). This means that for the field inside plasma on the plasma 
boundary the following condition is satisfied: 

e(-ai) = 65, e(+ai) = e^, (2.7) 

where ai = a/Hs the dimensionless width (depth) of slab. Further we will 
write a instead a\. 

We consider specular boundary conditions. For the function /i(x, \i) this 
conditions will be written in the following form: 

h{—a,iJ,) = h{—a,—iJ,), < /i < 1, (2.8) 

and 

h{a, iJ.) = h{a, —11), — 1 < /i < 0. (2.9) 

Condition of symmetry of boundary conditions (2.7) - (2.9) and the 
equation (2.6) mean, that electric field e{x) and function h{x,ii) in the 
slab possess properties of symmetry 

e{x) = e{—x), /i(x, /i) = /i(x, — /i). (2-10) 
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The non-flowing condition for the particle (electric current) flow through 
the plasma boundary means that 

1 1 

J iih{—a^ii)dii = J iih{a^ (i) dn = Q. 
-1 -1 
In accordace to (2.10) this condition it is carried out automatically. 
The problem statement is completed. Now the problem consists in flnd- 
ing of such solution of the system of equations (2.5) and (2.6), which sat- 
isfies the boundary conditions (2.7) and (2.8). Further, with the use of the 
solution of the problem, it is required to built the profiles of the distribu- 
tion function of the electrons moving to the plasma surface, and profile of 
the electric field. 

5. SEPARATION OF VARIABLES AND CHARACTERISTIC 

SYSTEM 

Application of the general Fourier method of the separation of variables 
in several steps results in the following substitution [31]: 

hr^{x, ii) = exp(-— )$i(77, ii) + exp(— )$2(77, lA^ (3-1) 

rj r] 

e,{x) = \eM-—) + eM—)]m. (3-2) 
L r] r] i 

where r] is the spectrum parameter or the parameter of separation, which 
is complex in general. 

We substitute the equalities (3.1) and (3.2) into the equations (2.5) and 
(2.6). We obtain the following characteristic system of equations: 

1 

zo{ri - i2)<^iiTj, ii) = riiJ,E{r]) + ^ J $1(77, /i')V, (3.3) 

-1 
1 

zoirj + Ai)$i(7], //) = wEir]) + | ^ $2(77, l^')dn', (3.4) 

-1 

1 

-'-^Eir^) = -,-l [ Mv.f^W. (3.5) 
rj 2 J 

-1 
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1 

rj 2 J 

-1 



Let us introduce the designations: 



2 ^ Zq 



V 



From equations (3.5) and (3.6) we obtain 



1 

/ 

-1 



$i(77,//)c?// + $2(^,/^) d^i = Q. 



Let us introduce the designations 



n{ri) = J ^i{r], iijdji. 



From equation (3.5) we obtain that 



whence 



or 



n{rf) = 



3 r] 

FJr, 

n{r]) = —2r]i 



n 



(3.6) 



(3.7) 



(3.8) 



(3.9) 



By means of equahties (3.7) - (3.9) we will rewritten the equations (3.3) 
and (3.4) 

(3.10) 



(77 - //)$i(77, 1^) = ^^il^r] - r]l), 



(3.11) 



Solution of the system (3.10) and (3.11) depends essentially on the con- 
dition whether the spectrum parameter rj belongs to the interval — 1 < 
77 < 1. In connection with this the interval — 1 < 77 < 1 we will call as 
continuous spectrum of the characteristic system. 
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Let the parameter 77 G (—1,1). Then from the equations (3.10) and 
(3.11) in the class of general functions we will find eigenfunction corre- 
sponding to the continuous spectrum 

Ml, = ^Pi^^^l^ + g,irj)6ir^ - ^), (3.12) 



r]-fi 



$2(77,/.) = ^pEl±li + ^,(^)^(^ _ (3.13) 

In these equations (3.12) and (3.13) 5{x) is the delta-function of Dirac, 
the symbol Px~^ means the principal value of the integral under integrating 
of the expression x~^. 

Substituting now (3.12) and (3.13) in the equations (3.5) and (3.6), We 
receive the equations 

1 



1-^ 
2c 



li'j] 



2 

m 



-1 



zq 



from which we obtain 



1 + ^ 
2c 



-1 



Ac 



77 



E{ri). (3.14) 



Here dispersive function is entered 



\{z) = 1 



z 

2c 



flZ 



2 

m 



(i/i. 



(3.15) 



where 



e Zq 



Functions (3.12) and (3.13) are called eigen functions of the continuous 
spectrum, since the spectrum parameter 77 fills out the continuum (—1, +1) 
compactly. The eigen solutions of the given problem can be found from 
the equalities (3.1) and (3.2). 
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The dispersion function X{z) we express in the terms of the Case dis- 
persion function [31]: 



1 1 z 



X{z) = l-- + -[l--)X,{z), 



where 



1 1 

dr 1 f T dr 



T — Z 2 J T — Z 
-1 -1 



is the Case dispersion function [31]. 

The boundary values of the dispersion function from above and below 
the cut (interval (—1, 1)) we define in the following way: 



The boundary values of the dispersion function from above and below 
the cut are calculated according to the Sokhotzky formulas 



A±(/.) = A(/.) ± ^{r^l - -1< < 1, 

ZTJi Zq 



from where 



ViZo 

^ = A(/i), -l<ii<l, 



where 



1 

-1 

and the integral in this equality is understood as singular in terms of the 
principal value by Cauchy. Besides that, the function A(//) can be repre- 
sented in the following form: 
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Substituting relations (3.14) in (3.12) and (3.13), we will present last 
expressions in the following form 



and 



or 



where 



^i(??,y") = 



^2(/7,/i) = 



P— — — - 2c^^6(t] - ii) 

Tj — jl 7] 



$1(77, /i) = ^F(/?,A^), $2(/?,/i) = ^F{-r^,ii)., 



and 



F(,.^) = P^-2o^i(,-p). 

Tj — jl T] 

It will be necessary for us the following relation of symmetry 

F{-r], -11) = -F{r],/i). 

Let us notice, that eigen function F(r], /i) satisfies to following condition 
of normalization 



J F{T],ii)dii 



-1 



= -2cP- 



and that analogous 



i 

/ F(-T],ii)dii = 2cP-. 
J 1 



So, eigen function of a continuous spectrum is constructed and it is 
defined by equality 

E{r]) 



hrjix^ji) = exp 
or, in explicit form 

/i^(x, //) = exp 



lF(7y,//) +exp (—]F[-r],ii) 



T] 



T] 



(- 



xzo \ iir] - r]{ 
f] J f] — fi 



+ exp 



xzQ \ idr] + r]f 
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-2c 



exp 



^^^]6{ri - /J,) - exp (^)S{ri + fi) 



Let us replace exponents by hyperbolic functions and also we will transform 
both square brackets from the previous expression. As a result we receive, 
that 



+2c^{-6{r] - ii) + 6{7] + ii)) +sh 
r] 



r]-ii 



r] 



+2c^{dir] - fi) + Sir] + fi)) 
We will designate further 

77 + /i 



7/ - /i 



= 2P 



rf — ji^ 



and 



V^(77, ]i) = -F{r], 11) + F{-r], ^z) = p( - 



T]- fi 



= 2P 



rf — 11^ 



Thus, eigen function of a continuous spectrum it is possible to present 
in the form of a linear combination of a hyperbolic sine and kosine 

u(r/2 - nl) X{r]) 



hn{x,/j,) = — ch 



P- 



rf — ji^ 



— c- 



r] 



— sh 



XZq 



jf — jl^ 



7] 



Let us notice, that eigen functions of a continuous spectrum it is possible 
to present and in such form 



hrj{x,ll) = 



E{r,) 



ch^(F(77,/i) + F(-7y,/i)) + 



14 



+ sh^(-F(77,//) + F(-7y,//)) 



(3.16) 



7] 

6. EIGENFUNCTIONS OF THE DISCRETE SPECTRUM 

According to the definition, the discrete spectrum of the characteristic 
equation is a set of zeroes of the dispersion equation 

^ = 0. (4.1) 

We start to search zeroes of the equation (4.1). Let us take Laurent 
series of the dispersion function: 

A(^) = Aoo + ^ + ^ + --- , \z\>l. (4.2) 

Here 

Aoo = A(oo) = 1 - — + ^ , 



A2 = 



zo SzqtiI 

1 /I 1 



zo\3 brjl'' 
1 /I 1 



zqVS 7r]l 

We express these parameters through the parameters Q and e: 

2{n -l)+ie + {n-l){n-l + ie) 



Aoo = A(oo) 



A, 

A4 = - 



9 + 5i£{Q + i£ 



i5{n + i£y ' 

15 + 7i£{Q + ie) 



36{n + iey 

It is easy seen that the dispersion function (3.15) in colhsional plasma 
(i.e. when £ > 0) in the infinity has the value which doesn't equal to zero: 

Aoo = A(oo) 7^ 0. 

Hence, the dispersion equation has infinity as a zero rji = oo, to which 
the discrete eigensolutions of the given system correspond: 

zq 
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This solution is naturally called as mode of Drude. It describes the 
volume conductivity of metal, considered by Drude (see, for example, [32]). 

Let us consider the question of the plasma mode existence in details. We 
find finite complex zeroes of the dispersion function. We use the principle 
of argument. We take the contour (see Fig. 1) 



which is passed in the positive direction and which bounds the biconnected 
domain Dr. This contour consists of the circumference 



and the contour 7£, which includes the cut [—1,+!], and stands at the 
distance of e from it. 

Let us note that the dispersion function has not any poles in the domain 
Dn. Then according to the principle of argument the number [33] of zeroes 
N in the domain Dr equals to: 



Considering the limit in this equality when £ — )• and taking into ac- 
count that the dispersion function is analytic in the neighbourhood of the 
infinity, we obtain that 




= : =R^ R = 



1 

-, e>0} 



N = 




1 



1 




-1 



-1 



1 




-1 



So, we obtained that 



1 




-1 




Fig 1. Contour = U 7^ for calculations of number of zero 
dispersion function. 
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We divide this integral into two integrals by segments [—1,0] and [0,1]. 
In the first integral by the segment [—1,0] we carry out replacement of 
variable r —r. Taking into account that A+(— r) = A~(r), we obtain 
that 

1 



If,, A+(r) If,, A+(r) 1 A+(r 
N = -— din — f4 = — T-fr = - arg ^ 



27rz J X (r) iri J A (r) tt A (r) 

1 



(4.3) 



A+(r) 

Here under symbol argG(r) = arg ^_ ^ ^ we understand the regular 

branch of the argument, fixed in zero with the condition: arg G{0) = 0. 
We consider the curve 

7 = {z: z = G{t), O^r^+1}, 

where 

It is obvious that 

G(0) = 1, lim G(r) = 1. 

T— >-+l 

Consequently, according to (4.3), the number of values N equals to 
doubled number of turns of the curve 7 around the point of origin, i.e. 

Ar = 2x(G), 

where 

x(G) = Ind[o,+i]G(r) 

is the index of the function G{t). 

Thus, the number of zeroes of the dispersion function, which are situated 
in complex plane outside of the segment [—1,1] of the real axis, equals 
to doubled index of the function G{t), calculated on the "semi-segment" 

[o,+i]. 

Let us single real and imaginary parts of the function G{ii) out. At first, 
we represent the function in the form: 



Giii) = 



{zo - l)r]l + {r]l - //2)Ao(//) - is{iJ,){r]1 - 
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where 



TT 



and 



is the dispersion function of Case, calculated on the cut (i.e., in the interval 

(-i.i))- 

Taking into account that 

-2 .en 



Zq — 1 = —I— = — Z — , 

v e 



2 

^1 = ^ 



.en 



we obtain 



where 



{zo - 1)7^1 = -— - i— 
p-{li)+iQ-{li) 



P+ili)+iQ+ili) 



Q^i/l) = £n{l + Ao(yw)) ± S{ll){£^ - 3//2). 

Now we can easily single real and imaginary parts of the function G{fi) 



out: 



9W 9KN 



Here 



g{li) = [P^{ii)f + [Q^{ii)f = [n" + Ao(3a.^ - e')- 
-ensf + [en{l + Ao) - 5(3^^2 - e^)f, 

-e'n'is' - (1 + Xof] - (3Ai' - sYs^ 
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We consider (see Fig. 2) the curve L, which is defined in imphcit form 
by the following parametric equations: 

L = {{Q,e): gi{iJ,;Q,e) = 0, p2(/^; = 0, < // < 1}, 

and which lays in the plane of the parameters of the problem {^,e), and 
when passing through this curve the index of the function G{ii) at the 
positive "semi-segment" [0, 1] changes stepwise. 
From the equation g2 = we find: 

Now from the equation gi = with the help of (4.4) we find that 

e = /L^, (5.5) 
where „ „ 



AoM[«2(p) + (l + Ao(p))2] 
Substituting (4.5) into (4.4), we obtain: 



n = ^^jLiiii), (4.6) 
3/i2[s2(^) + Ao(^)(l + Ao(/i))]2 



where 

Ao(/.)[s2(/.) + (l + Ao(/.))2] • 
Functions (4.5) and (4.6) determine the curve L which is the border if 
the domain (we designate the external area to the domain as D~) in 
explicit parametrical form (see Fig. 2). As in the work [34] we can prove 
that if (7,5) e ^ then niG) = Ind[o^+i]G'(/i) = 1 (the curve L encircles 
the point of origin once), and if (7, e) G then x(G) = Ind[o,+i]G(//) = 
(the curve L doesn't encircle the point of origin). 

We note, that in the work [34] the method of analysis of boundary 
regime when (n,£) G L was developed. 

From the expression (3.2) one can see that the number of zeroes of the 
dispersion function equals to two if {^^e) G ^ and equals to zero if 
iSl.e) e D-. 
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Fig. 2. 



Since the dispersion function is even its zeroes differ from each other by 
sign. We designate these zeroes as following itr/o, by 770 we take the zero 
which satisfies the condition Re 770 > 0. The following solution corresponds 
to the zero rjQ 

XZq 

hr,o{x, 11) = ch [F{r]o, //) + F{-r]o, //)]+ 

^0 

+ sh^[-F(r/o,Ai) + F(-^o,//)], (4.7) 
^0 

e,,(.;) = 2ch^. (4.8) 

Here 

r]o- Id tlQ^ 11 

It is easy to see, that function hrj^i^x, ji) is even on 770 
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Function /i,7o(x, fi) we will present in the explicit form 



hrjo{x,ii) = ch 



XZq 



+ 



+ sh 



XZq 



or 



hrjoix, fi) = ch ^(p{t]q, ii) + sh ^V^(77o, m)- 

^0 ^0 



Here 



(4.9) 



(4.90 



^'(^o, IJ) = -F{r]o, 11) + F(-7yo, /^) = 



2a^(^o - vi) 
9 9 
^0 



This solution is naturally called as mode of Debay (this is plasma mode). 
In the case of low frequencies it describes well-known screening of Debay 
[3]. The external field penetrates into plasma on the depth of vd, 
is the raduis of Debay. When the external field frequencies are close to 
Langmuir frequencies, the mode of Debay describes plasma oscillations 
(see, for instance, [3, 32]). 

Note 5.1. If to enter expression c/ z "inside" of expression of dispersion 
function we will receive expression for dispersion function h{z) from our 
article [21] 

Kz) = '-X{z) = '--^-]^^d^ = 
^ ^ z ^ ' z 2j fi-z ^ 

-1 

= z — {z — 777)- m . 

z ^ '^'2 z + l 

6. EXPANSIONS IN THE TERMS OF EIGEN FUNCTIONS 
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We will seek for the solution of the system of equations (2.5) and (2.6) 
with boundary conditions (2.7)-(2.9) in the form of linear combination 

of discrete eigen solutions of the characteristic system and integral taken 
over continuous spectrum of the system. Let us prove that the following 
theorem is true. 

Theorem 6.1. System of equations (2.5) and (2.6) with boundary con- 
ditions (3.1), (3.6) and (2.7) has a unique solution, which can be presented 

as an expansion by eigen functions of the characteristic system: 



ch —{F(n^, ii) + F(-77o, //)) + sh — (-F(7yo, IJ) + F(-ryo, //)) 



+ 



1 



+ 




ch 



-1 



7] 



F{r],ii)+Fi-Ti,ii) 



+sh 



XZq 



-F(77,//)+F(-77,Ai)l j^rfT], 



1 



ZqX 



e(x) = Eoo + 2Eoch ^ + 2 

r] 



ch -^—Eirj) drj. 



(6.1) 
(6.2) 



Here Eq and Eoo is unknown coefficients corresponding to the discrete 
spectrum {Eq is the amplitude of Debay, Ei is the amplitude of Drude), 
E{rj) is unknown function, which is called as coefficient of discrete spect- 
rum. 

When (0,£) G D in expansions (6.1) and (6.2) we should take Eq = 0. 

Further we will consider the following case (Jl, e) G Z)+. 

Our purpose is to find the coefficient of the continuous spectrum, coeffi- 
cients of the discrete spectrum and to built expressions for electron distri- 
bution function at the plasma surface and electric field. 

Proof. Let us consider expansion (6.1), we will replace in it ^ on — /i. 

1 

Then we will substitute h{a, fi) and h{a, —fi) in The equation -[hia^fi) — 
h{a^ — //)] = 0. After variety of transformations let us have 



E^fi + Eo ch — F(7yo, + F{-r]Q, ji) 



+ 
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+ / ch ^[F(ry, fi) + F(-77, fi)]E{r])dr] = 0, -!<//<!, (6.3) 
J r] 



-1 

or 



Eoofi + £^0 ch ^(/p(t/o, /i) + 2 / ch ^F(7y, ii)E{r])dr] = 0. (6.4) 

-1 

Let us pass from integral Fredholm equation (6.4) to singular integral 
equation with Cauchy kernel, having substituted in (6.4) obvious represen- 
tation F{r], ji) 

1 

E^li + E,ch ^^(^o,A^) + 2 / ^^!^^(^)ch ('-^)dr^- 

-1 

-4c^E{fi) ch ^ = 0, 0<fi<l. (6.5) 

It is easy to check up, that function 

1 

M{z) = / ^^^^E(/7) ch ^df] (6.6) 

J f]- Z 7] 

-1 

is odd. Besides, all members of the equation (6.5) are odd on fi. 

Extending the function E{r]) into the interval (—1,0) evenly, so that 
E{r]) = E{—r]), and extending the equation into the interval (—1; 1) un- 
evenly, we transform the equation (6.5) to the form 

1 

EooA^ + EoCh ^^(770,A^) + 2 / Ch (^)^^^E(77)d77- 

-1 



4c^ch (^)e(;^) = 0, -1<a.<1. (6.7) 



Let us reduce the equation (6.7) to boundary value problem of Riemann 
— Hilbert. For this purpose we will take formulas of Sohotsky for the 
auxiliary functions M{z) and dispersion function X{z) 

M+{fi) - M-{fi) = 27nch (^^Yfi'^ -r]l)E{fi), -1 < fi< 1, (6.8) 
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M+(//) + M-(//) 



-I < ll<l, 



where 



and last integral is understood as singular in sense integral principal value 
on Cauchy, and 



ITT 



A+(/.)-A-(A.) = -^M/^^-^f): 



A+(//) + A-(//) 



-1 < // < 1, 

-1 < yLi < 1, 



where 



-1 

As a result of use of last formulas we will come to the boundary value 
problem 

Eoofi + Eoch ^^(77o,Ai)+ [M+{fi) + M-{fi)] + 



A+(/i) + A-(//) 



[M+ifi) - M-(//)] = 0. 



A+(^)-A-(Ai) 
Let us transform this equation to the form 

(M+ + M-){X+ - A") + (M+ - M-)(A+ + A-)+ 



azo 



^0 



-1 < li<l. 



From here we receive the boundary condition of boundary value problem 
of Riemann — Hilbert 

A^-(;.) [m+(/.) + %/i + f eh ^^(ryo, f^) 
L 2 2 770 

= A-(A.)[M-(;L.) + %A^ + ^ch^v:.(/7o,/i)] = 0, -1< /i < 1. (6.9) 
Let us copy this problem in the form 



- = 0, -1< // < L 



(6.10) 
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In the problem (6.10) $^(//) is the boundary values in interval — 1 < // < 1 
of function 

r 1 1 

^z) = X+{z) M+{z) + -E^z + -Eo ch ^(/9(77o, z) 

which is analytical function in complex plane with cut C \ [— 1, 1]. 

The problem (6.10) is the problem special case about jump. It is prob- 
lem of definition of analytical function on its jump on a contour L: 

$+(//)-$+(//) = ^(//), fieL. 

The solution of such problems in a class decreasing at infinitely remote 
point of functions is given by integral of Cauchy type 

1 f (p{T)dr 



^z) = 



27Ti J T — Z 
L 



However, in the problem (6.10) unknown function ^{z) has at infinitely 
remote point following asymptotic 

^{z) = 0{z), z ^ oo. 

Therefore (6.11) it is necessary to search for the solution of the problem 
in the class of the growing as z at vicinity of infinitely remote point. 

According to [33] the general solution of the problem (6.10) is given by 
the formula $(z) = Ciz. 

In explicit form the general solution of the problem (6.10) registers as 
follows: 

r 1 n.zn 1 1 



X{z) M{z) + i^och —ifim, z) + ]-E^z 
z r]Q z . 



where C\ is the arbitrary constant. 

From this general solution we can find function M{z) 

Miz) = -\e^z - \e, ch ^v^(77o, z) + ^ (6.11) 

Let us remove the pole at the solution (6.11) at infinitely remote point. 
We receive, that 

Ci = \e^\oo. (6.12) 
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8. COEFFICIENTS OF THE CONTINUOUS AND 
DISCRETE SPECTRA 

Now we will eliminate polar singularity at the solutiom (6.11) at points 
±77o- Let us allocate in the right part of the solution (6.11) members 
containing the polar singularity at points z = 770. In the point vicinity 
z = rjo taking into account equality A (770) = fairly following expansion 



Miz) = -lEo,z-lEoch^ 

Cm 



+ 



From here is visible, that for pole elimination at the point z = 770 is 
necessary to equate to zero expression in a square bracket, calculated at 
z = rjQ. Then we receive, that: 

m ° A'(77o)(77? -/??) 

The coefficient of continuous spectrum E(r]) we will find from the for- 
mula (6.9) 

azo 1 M+{r]) - M-{r]) 



ch ^Eq = -^jYzrxy^ 2V- (7-1) 



E{ti) ch - 2 
77 ziri Tj'^ — rii 

^ooAoo J] f I 1 



2 27rz(r/2 - 7/f) VA+(r/) \-{t]) 
Let us take advantage of the boundary condition for electric field 

1 



(7.2) 



Eoo + 2Eo ch — + 2 / Eiri) ch —dr] = e,. (7.3) 

Vo J V 
-1 



Let us substitute coefficient of the continuous spectrum (7.2) and coef- 
ficient of discrete spectrum (7.1) in the equation on electric field (7.3). We 
will receive the following equation 



1 



-1 
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or 



E A 



2m 



rjdr] 

Aoo >^'iv)iVo-vl) ' 27TiJ yX+ir]) \-{iri)) rf -7)1 

-1 



+ 



Integral from this expression we will calculate by means of the contour 
integration and the theory of residual 

1 




1 



1 



A+(77) A (77)/ rf — 



ijdr] 



Res + Res + Res + Res + Res 

L 00 -771 rji rjo -??o 



Z 



Kz){z^-Vl) 



1 1 

+ ^ + 



2r/o 



Aoo ■ Ai ■ x'm^-r^ir 

Substituting this integral in the previous equation, we find that 

Ai 



A 



(7.4) 



00 



Now by means of a relation (7.1) we find 



Eoch 



-e. 



(7.5) 



m 'X'{iio){r]^-r]l)' 

By means of the found coefficients of discrete and continuous spectra 
we find an electric field profile in the slab 

e{x) Ai 2r/oAi 



Ai 



27ri 



-1 



Aoo A'(77o)(77o -^1) 

1 1 1 ch [x zq / ij) ijdr] 

A+(77) X~ {j]) } ch [azo / T]) T]'^ — j]i 



(7.6) 



Integral from (7.6) we will calculate by means of residual and contour 
integration 

1 



1 



27ri 



1 



1 



-1 



ch(xzQ/r]) T]dT] 
X~^{r]) X~ [t]) 1 ch [azQ / T]) T]"^ — j]l 
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+ 00 

Res +Res +Res +Res +Res + \^ Res 

00 -rji rji rjo -rjo *fc 

fc= — GO 



2;ch [xz^j z) 



sh {az{)l z){z^ — tjDX^z) 



1 ^ ch{xzo/r]i) ^ 2r]och{xzo/r]o) 



Aoo Ai ch (azo/ryi) A'(7yo) W " Vi) ch {azo/r]o) 



E 



Here i]^ are zeros of function ch(azo/2;), 

also it has been considered that sh(a2;o/tfc) = i(— 1)^. 

Hence, an electric field profile in the slab of degenerate plasma is possible 
to express without quadratures 



e{x) _ ch (3:2:0/7^1) _ ^ A (^W^fc) 

ch (azo/771) az^ Ktk){tl - ril) sh (azo/tk) ' 



A; 

or, that all the same, 



ejx) _ ch {xzo/t]i) ^ ^ ch (xzp/tfc) 

eg ch(a2;o/77i) a^o , Ktk){tl-T]l) 



X 



Here ch [xz^/tk] = ch — (tt/c + 7r/2). 

a 

Integrating expression (7.7) on x from —a for a, we receive 
Q,= I e{x)dx = e,^ th ^ - e,^ ^ —| ^. (7.8) 

Thus, all the coefficients of the expansions (6.1) and (6.2) are found 
single- valued: the coefficient £"00 is determined according to (7.4), the co- 
efficient Eq is determined according to (7.5), the coefficient of the con- 
tinuous spectrum E{r]) is determined according to (7.2). Finding of the 
coefficients of discrete and continuous spectra of the expansions (6.1) and 
(6.2) completes the proving of existence of this expansions. Uniqueness of 
the solution in the form of expansions (6.1) and (6.2) can be proved easily 
with the use of contraposition method. 
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9. ENERGY ABSORPTION IN THE SLAB 

The energy of an electromagnetic wave absorbed in the slab of degen- 
erate plasma is calculated under the formula 

a 

Q = -^ [ lmE{x')dx'. (8.1) 

OTT J 

—a 

Here dimensional quantities E{x),x,a are written, and E(0) = Eq. 
Let us pass in (8.1) to dimensionless quantities 

It is as a result received 

a /I 

Q = - I t \2 ^^ / e(xi)dxi. (8.2) 

-a/l 

Further we will designate 

a /I 

Qi= J e{xi)dxi. 

-a/l 

Then according to (8.2) we have 

Q = ^^-^^lr.Q. (8.2') 

Besides, quantities xi and ai = a/l we will designate again through x 
and a. 

Earlier for electric field which we will present in the form 

1 

e{x) = + 2Eo ch ^ + / ch —E{r])dr], (8.3) 

-1 

coefficients of continuous and discrete spectra have been calculated in the 
explicit form. 
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Integrating (8.3), we receive 

1 

Qi = 2aEoo + 4Eo^ sh ^ + - / rjEirj) sh ^ drj. (8.4) 

-1 

Substituting in (8.4) coefficients of discrete and continuous spectra, we 
receive the following expression for electric field 

e{x) Ai 4Ai7yoth(a2;o/r^o) , 



2Ai 1 
zq 27ri 



-1 



rf th {azQ/ri) 
rf — rji 



(8.5) 



Integral from last expression we will calculate by means of the contour 
integration and the theory of residual 

1 



1 



27ri 



1 



1 



-1 



A+(77) A+(77) 



Tf ill {(IZo/t]) 



rf — rji 



Res +Res +Res +Res +Res + \^ Res 

00 -ryi 771 -770 r]o ^-^ tu 



k=—oo 



th {az{)l z) 



+00 



/4 



^^!h^Yi^+ ^^0 (q^o/^o) L 

^00 Ai 771 X'{r]o){r]^-r]l) azo X{tk)itl - r]l)' 



(8.6) 



where points tk are entered earlier. 
Substituting (8.6) in (8.5), we receive 

Qi 2771 a2;o 2Ai 
— = th - 



-00 



az^ ^ 



k=—oo 



(8.7) 



that is equivalent to expression (7.8). 



10. CONCLUSION 



In the present paper the linearized problem of plasma oscillations in slab 
(particularly, thin films) in external longitudinal alternating electric field 
is solved analytically. Specular boundary conditions of electron reflection 
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from the plasma boundary are considered. Coefficients of continuous and 
discrete spectra of the problem are found, and electron distribution func- 
tion on the plasma boundary and electric field are expressed in explicit 
form. 

Separation of variables leads to characteristic system of the equations. 
Its solution in space of the generalized functions allows to find the eigen 
solutions of initial system of equations of Boltzmann — Vlasov — Maxwell, 
correspond to continuous spectrum. 

Then the discrete spectrum of this problem consisting of zero of the 
dispersion equations is investigated. Such zeros are three. One zero is infi- 
nite remote point of complex plane. It is correspond to the eigen solution 
"Drude mode" independent of boundary conditions. Others two points of 
discrete spectrum are differing with signs two zero of dispersion function. 
These zero are correspond to the eigen solution named " Debay mode" . 

It is found out, that on plane of parametres of a problem {^,e), where 
Q, = uj/ujp, e = ly/ujp, there is such domain (its exterior is domain D~), 
such, that if a point (Q, e) G D~ Debay mode is absent. 

Under eigen solution of initial system its general solution isobtained. 
By means of boundary conditions in an explicit form expressions for co- 
efficients of discrete and continuous spectra are found. Then in explicit 
form (without quadratures) absorption quantity of energy of electric field 
in slab of degenerate plasmas is found. 
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